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We derive the classical transport equation, in scalar field theory with a 
g 2 V((j)) interaction, from the equation of motion for the quantum field. We 
obtain a very simple, but iterative, expression for the effective action T which 
generates all the ra-point Green functions in the high-temperature limit. An 
explicit and closed form is given for T in the static case. 



I. INTRODUCTION 

In thermal field theory, one is often interested in the behavior of Green functions at high 
temperatures. These hard thermal functions have been much studied in QCD ]T|-0 where 
they have a fundamental role in the resummation procedure of Braaten and Pisarski pi. 
This program is usually based on quantum field theory, which involves in general rather 
complicated calculations. On the other hand, it was shown that the hard thermal effects 
can be described in a simpler and more transparent way using classical transport equations 

One may ask about the relevance of considering this approach in the framework of a scalar 
field theory, given that it has been already applied to more realistic theories like QED and 
QCD. Our motivation is twofold. First, we recall that an essential ingredient in the derivation 
of classical transport equations is the equation of motion for a particle propagating in a 
external background field. In QED, the equation of motion for the electron is governed by 
the well known Lorentz force. The natural generalization to a non-Abelian gauge theory like 
QCD, is the Wong equation [ 10| . The geodesic equation in gravity ]TlJ] is another example 



where the particle has a well defined equation of motion. The common feature of these 
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theories is that they are based on a gauge invariance principle which prescribes in a well 
defined way how the particle couples to the external field. Nonlinear scalar theories, on the 
other hand, have no such a simple prescription and require a rather different approach [for 



an alternative treatment, see also reference [12]]. Our second purpose is to obtain an exact 



non-perturbative expression for the effective thermal action, which contains the dominant 
static contributions arising in each order of perturbation theory. (Such a closed form result 
for the corresponding action in gauge theories is generally not known.) 

In section 2 we derive the Boltzmann equation, by considering the expectation values of 
the relevant quantum-mechanical quantities. The analysis shows that the condition for the 
classical limit requires that the distances over which the background field varies must be 
large compared with the particle de Broglie wavelength, i.e: 

1 ft / 1N 

(l) 

where A; is a typical wavenumber of the background field and p is the kinetic momentum of 
the particle. 

The behavior of Green functions at high temperatures is governed by the region involving 
large values of momenta of the particle propagating in the thermal loops, a domain which is 
consistent with the condition given by Eq. ([!]). Using the properties of the classical transport 
equation, we show in section 3 that the effective action which generates the high-temperature 
contributions of the Green functions in this domain has a form such that: 

d ^ 



5V( 2 )((j>) 



l = 



2p ■ d ^ ^ ^ dp< 



Hp) , (2) 



where Ct n ) is an appropriate normalization constant which depends on the space-time di- 
mension n, and fo(p) is a function proportional to the Bose-Einstein probability distribution. 
This result is not obvious, so far as we know, since in scalar field theory the n-point Green 
functions are sub- leading for n > 2. 

The above action has two interesting features: 

• a non-locality in configuration space in the form of powers of the operator (p ■ d)^ 1 . 



• in the static case the effective action can be expressed in closed form in terms of 
modified Bessel functions. 



II. DERIVATION OF THE BOLTZMANN EQUATION 



Consider the scalar field theory described by the Lagrangian density: 

where m is the mass of the field and V(<E>) describes the interactions between the scalar 
fields. Examples of particular interest are the $ 4 theory in 4 space-time dimensions and the 
$ 3 theory in 6 dimensions. In what follows we shall work in units where c = 1, but will keep 
explicit the dependence on the Planck constant. 

We now divide $ into a sum of a quantum field ip, associated with the particle, and a 
classical background field (f>: 

$ = if; + . (4) 
For our purpose, it is sufficient to consider only the terms in C which are quadratic in ip: 

c2 = -\ o^r - ^ ^ - ^ . (5) 

This leads to the linearized Euler-Lagrange equations: 

h 2 d 2 i)= [m 2 + g 2 V {2 \(f))}^. (6) 

It is convenient to write this equation in Hamiltonian form. To this end we introduce 
two wave functions ipi and tp2 [ 13f 



dip 

i) = + ip 2 , mtyi-fo) =ih — , (7) 
together with the single-column matrix \1/ and the Pauli matrices given by: 
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It is then straightforward to verify that (^) is completely equivalent to the equation of 
motion: 



ih 



~dt 
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where the Hamiltonian operator H is given by: 

H = ±(cr 3 + ia 2 )[P 2 + g 2 V^ 2 \<P)} + mcr 3 , 
2m 

and P is the canonical momentum operator —ihV. 

In the Heisenberg representation, we have the operator equation: 



(ii) 
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Using the form of the Hamiltonian given by (JTTJ), this leads to the equation of motion 
for the canonical momentum: 



dP 

~dt 
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+ ia 2 )S7V^\<P). 
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We now consider the expectation value of this quantity: 

dP r „ n „* dP 



< — > = / d n l X ^T(T3 — * 

dt J dt 



(14) 



To evaluate the space integral, we make the assumption that the background field is 
roughly constant over a distance of the order of the particle de Broglie wavelength. Taking 
into account the spatial Lorentz contraction, we then obtain: 



(15) 



where r is the proper time of the particle. 

Let us compare (15) with the classical equation of motion of a scalar particle moving in 
an external field, which is obtained varying the action: 
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Proceeding as usual, we get the following equation of motion: 
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dr 
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where p^ denotes the kinetic momentum of the particle. Note that ( |T7D is consistent with 
the on-shell condition p 2 = —m 2 which should hold for a physical particle. 

In order to compare this relation with equation fllBD , we replace p ■ d by m dV^ / dr, 
which follows by the use of the chain rule. Then, we can write ( |17|) in the form: 
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This expression is in agreement with fllBD , provided the canonical and kinetic momentum 
are related by: 
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Furthermore, the classical limit requires the condition: 



(20) 



where A; is a typical wavenumber of the background field. 

Equation ([[7]) can be used to determine the evolution of the distribution function f(x, p), 
which represents the probability of finding a particle in the state (x, p). From the property: 

df(x, p) 



dr 







(21) 



which holds in the collisionless case, and using the constitutive relation 
p^ = mdx^/dr, we obtain the Boltzmann equation: 



P -9f(x, p) - 



2 l-g 2 V^/2p 2 



p* 
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dp 11 
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We remark that the Planck constant does not appear in this equation, as one would have 
expected for a classical transport equation. 

Let us expand f(x, p) in powers of the coupling constant: 



f(z, P) = fo(p) + 9 2 fi(x, p) + g 4 f 2 (x, p) + • • 



(23) 



where: 



fo(p)=0(Po)5(p 2 + m 2 )N(p o ). 



(24) 



Here, the theta and delta functions guarantee positivity of the energy and on-shell evo 
lution, and N(p ) is the Bose probability distribution: 
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Inserting the expansion of f(x, p) into the transport equation (|22|) and identifying the 
coefficients of the powers of the coupling constant, we can determine recursively the com- 
ponents fi(x, p). We find, for example: 



fi(x, P) = 2 



p ■ d 



P 



dfo 



(26) 
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where • • • denote similar terms involving powers of 1/p 2 . It turns out that such terms are not 
connected with the contributions obtained, in the high-temperature limit, from the 1-particle 
irreducible Green functions. 
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III. THE EFFECTIVE ACTION 



In order to determine the form of the effective action, we will consider only those terms 
in the expansion of f(x, p) which involve powers of the operator: 

which has in configuration space a non-locality of the form (p-d)^ 1 . Note that in momentum 
space, D is a function of degree zero in k. 

The above terms are generated by the series: 

oo 

F{x,p) = Y,9 2l D l {4>)h(p)> (29) 
1 = 

which satisfies the differential equation: 

p.dF(x,p)-£ (0) = . (30) 

It is next shown that this function is relevant in connection with the generating functional 
T defined by: 

oo „ 

= C [n) g 2 Y.9 21 ctpD l (<i>)f (p), (31) 
i = o J 

where Cr n ) is a normalization factor. 

To illustrate the content of the above quantity, let us consider for definiteness the <f) A 
theory in 4 space-time dimensions. For I = 0, we then get in momentum space the expression: 

which represents the scalar self-energy function II2 shown in Fig. 1, provided we choose 
C(4) = l/(27r) 3 . As is well known [14] this function has a leading behavior proportional to 



rp2 

For I = 1, a simple calculation which makes use of integration by parts, gives in momen- 
tum space a contribution like: 



sm = -ik? I d,p9[Pa)S[p2 + m2)NM ww • (33) 

where Tik denotes the total momentum of a pair of fields 2 . 

We shall now compare this result with the one obtained from the 1-particle irreducible 
Green function shown in Fig. 2. This contribution can be expressed in terms of a forward 
scattering amplitude of an on-shell particle in an external field || as indicated by the two 
graphs in Fig. 3. Also, an integration over the particle's momenta p^ must be performed 
with a weight factor given by the Bose probability distribution N(p ). We then obtain a 
contribution of the form: 

1 
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(34) 



_h 2 k 2 + 2hk-p Ti 2 k 2 -2Tik-p 
We now consider the high-temperature behavior, which is governed by the region involv- 
ing high values of p. Expanding the Feynman propagators in this region, we get: 

1 1 k ,~ , Ti k , , , 

h 2 k 2 + 2hk- P + h 2 k 2 -2hk- P = ~2( P -k) 2 + ^~p^'~ ( ' 

Note that the leading term, which is of degree zero in k, is classical since the % dependence 
cancels out. Therefore, as far as the high-temperature domain is concerned, (0) effectively 
reduces to the classical amplitude (|33D, yielding a InT contribution. 

These calculations show that the operator D generates iteratively, in the high- 
temperature limit, the dominant contributions associated with the 1-particle irreducible 
Green functions. It is clear that this mechanism generalizes to more general forms of inter- 
actions between the scalar fields. Hence, the generating functional T given by equation fl3~T|) 
can be interpreted as the high-temperature effective action in scalar field theory. 

We have not been able to evaluate T in closed form for arbitrary external fields. How- 
ever, this is possible when the fields are static, in which case the function F(x, p) can be 
determined exactly. It is easy to verify that the static solution of the equation fl30D , which 
reduces to fo in the absence of interactions, is given by: 

F(x, p) = 6( Po )5(p 2 + m 2 + ^V (2) (0))iV(p o ) . (36) 



With this form, (pTD may be integrated to give the following local expression: 

r = C (n) J d n x J d n P 6(p )6[-(p 2 + m 2 + <?V (2) (0))MPo) . (37) 

Using the constraints imposed by the 9 functions, the p integrations can be done, with 
the result: 

/ d n x r dp, (p 2 - M 2 Y n -^ 2 N(p ) , (38) 



2C [n) n^ 2 



n-1 r[(n - l)/2] 
where M 2 = m 2 + g 2 V {2) ((p). 

After expanding N(po) as a power series in exponentials, the po integrations can be 

expressed in terms of modified Bessel functions |15[]. We then obtain, in the static case, the 

following expression for T: 

r = ^C {n) f (^) W2 / rxM^K n/2 (^j . (39) 

The properties of this result may be illustrated in the case of the 4 theory in 4 space- 
time dimensions. Making a high-temperature expansion of the Bessel function K 2 in powers 



of IM/T, and dropping an irrelevant constant term, we obtain a series of the form: 



T 



2 



d 4 xM 2 



T + 8 T 2 " n T 2- + 24T 4 + 



(40) 



where £W i s the derivative of the Riemann zeta function and 
M 2 = m 2 {l + g 2 (f) 2 /2m 2 ). 

The first term in equation ([5]) corresponds to the leading T 2 contribution associated 
with the scalar self-energy function. The second term contains, in addition to sub-leading 
self-energy function, the InT contribution of the 4-point function in the high-temperature 
limit. The contributions associated with the scalar 6-point function first appear in the third 
term. 



It is interesting to remark that the expression [40] also contains a non-perturbative term 
involving ln[l + g 2 <p 2 /2m 2 }. This contribution may be expanded in a power series of the 
coupling constant when g 2 cf) 2 < 2m 2 , but in very strong external fields the perturbative 
expansion is no longer meaningful. Hence, aside from the static condition, we expect that 
the result ( |39| ) may be applicable in perturbation theory for not too strong external fields. 
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FIGURE CAPTIONS 



Fig. 1 The scalar self-energy function n 2 . 
Fig. 2 The irreducible 4-point function II4. 

Fig. 3 The forward scattering amplitudes associated with the 4-point function. 




p p + hk p p p — hk p 



FIG. 3. 
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